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Abstract 

We consider sequences of rational interpolants r„(z) of degree n to the exponential 
function associated to a triangular scheme of complex points {zj^"^}j!iQ, n > 0, such 

that, for all n, \z^^'^\ < cn}~'^, j = 0, . . . , 2n, with < a < 1 and c > 0. We prove 
the local uniform convergence of r„(z) to in the complex plane, as n tends to infinity, 
and show that the limit distributions of the conveniently scaled zeros and poles of r„ are 
identical to the corresponding distributions of the classical Pade approximants. This ex- 
tends previous results obtained in the case of bounded (or growing like log n) interpolation 
points. To derive our results, we use the Deift-Zhou steepest descent method for Riemann- 

(2n) 

Hilbert problems. For interpolation points of order n, satisfying \zj \ < cn, c > 0, the 
above results are false if c is large, e.g. c > 27r. In this connection, we display numerical 
experiments showing how the distributions of zeros and poles of the interpolants may be 
modified when considering different configurations of interpolation points with modulus of 
order n. 

AMS classification: 41A05, 41A21, 30E10, 30E25, 35Q15 
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1 Introduction and main results 

Rational approximants to the exponential function have been the object of numerous studies in 
the literature. One motivation comes from the fact that the approximation of the exponential 
function naturally appears in many problems from applied mathematics, like, for instance, the 
stability of numerical methods for solving differential equations, the modeling of time-delay 
systems to be found, e.g. in electrical or mechanical engineering, and the efficient computation 
of the exponential of a matrix. Another, more theoretical, motivation comes from the partic- 
ular properties of the exponential and its approximants in the framework of function theory. 
One classical example of such properties is Pade's theorem about the convergence of Pade ap- 
proximants to the exponential, and its connection with deep results in analytic number theory. 
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Another typical problem has been the one of finding the rate of rational approximation to the ex- 
ponential on the semi-axis, the so-called 1/9 -conjecture. It attracted the efforts of many authors 
in the eighties and was eventually solved in [fTO l. 

The behavior of Pade approximants to the exponential function has been studied, among 
others, in [[T4l[T5l[T6l[T9l l9l. and for extensions to Hermite-Pade approximants, one may consult 
inillSliniini. Generalizations to rational interpolants are investigated in ll3l l4l l2l IIUI IITI . In 
it is shown that rational interpolants to the exponential function with bounded complex 
interpolation points (also with points growing like a logarithm of the degree) converge locally 
uniformly in the complex plane, as the degree of the interpolant tends to infinity. The proof 
uses the Deift-Zhou steepest descent method for Riemann-Hilbert problems [|8l|51|6ll3. In the 
present paper, we consider the case of interpolation points whose modulus may grow with the 
degree n of the interpolants, namely like n^~", < a < 1, and we show that the Deift-Zhou 
method can still be used to show convergence of the interpolants as n — )■ oo. For interpolation 
points whose growth is linear with respect to the degree, it is easy to see from the periodicity 
of the exponential on the imaginary axis that convergence cannot always hold true. Also, for 
the particular case of shifted Pade approximants, interpolating the exponential at the point n^, 
^ G C, it is possible to give a necessary and sufficient condition on ^ for convergence to hold 
true, see fi22J . 

Let us now describe our findings in more detail. Given a triangular sequence of complex 
interpolation points rii + ns > 0, we are interested in the behavior, as rii + n2 

becomes large, of the rational function — , with p„j , g„2 polynomials satisfying the conditions: 



(i) degpni < ni, degg„2<n2, (1.1) 



(ii) en,,n2{z) ■■= Pm{z)e + qn2{z)e'/^ = C(w„,+„2+i(z)), (1.2) 

(711+712) ■ n 1 

as z ^ Zj ' , J = 0, . . . ,ni + n2, 



with 



711+712 

a;n.+„2+iW= n (^-#^^"^^ (1.3) 

j=0 

For any choice of (possibly multiple) interpolation points, nontrivial polynomials p„j and qn2, 
such that ( |l.l[ )-( 1-2) hold true, always exist since these conditions are equivalent to a system of 



2n + 1 homogeneous linear equations with 2n + 2 unknowns. 
In this paper, we will only be interested in the diagonal case 

degpn<n, degg„<n, (1.4) 

and 

Pn{z)e-'^^ + qn{z)e'/^ = O{c02n+iiz)), z^zf''\ j = 0,...,2n, (1.5) 

though the general case could be studied similarly. As we will see in the sequel, even if we 
restrict ourselves to the diagonal case, pairs of polynomials of type (n — 1, n + 1), that is of 
degrees respectively less than or equal to n — 1 and n + 1, will show up in the study. 
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Let us write 



p„ :=max{|2f"^| : J = 0, ...,2n}, neN. (1.6) 

If the interpolation points do not grow too rapidly with n, i.e. if there exist constants < a < 1 
and c > such that 

1 — a 

Pn< o log n + c, n G N, (1.7) 



it was proved in [21, Theorem 2.2] that a pair {pn,qn), such that (1.4)-(1.5) hold true (with 
111 = n2 = n), satisfies 

locally uniformly in C, where g„ is normalized so that g„(0) = 1. In particular, pn{z) / qn{z) 
converges to — uniformly on compact sets in the complex plane as n — ^ oo. Our aim is to 
weaken the condition to interpolation points for which there exists < a < 1 and c > 
(independent of n) such that 

Pn < cn^^'^, neN. (1.8) 

This is our main result. 

Theorem 1.1 Let z'f'^\ n > 0, j = 0, . . . , 2n, be a family of interpolation points satisfying 
( |i.<g| ) with < a < 1 and c > 0. Let pn and qn be polynomials satisfying {L4)-{L5). Then, the 
following three assertions hold true: 

(i) All the zeros and poles of Vn = Pn/qn tend to infinity, as n becomes large, and, more 
precisely, no zeros and poles of Tn lie in the disk {z,\z\ < pn}, for n large. In partic- 
ular, dividing equation {L5) by qn, we get r„ = Pn/qn <^ rational interpolant to — 
satisfying 

+ rn{z) = 0{u2n+i{z)), as z ^ zf''\ j = 0,. . . ,2n. 

( ii) As n ^ CO, 

p^{z)^-e''\ qn{z)^e-''\ ^{z) ^ -e% (1.9) 
locally uniformly in C, where qn is normalized so that g„(0) = 1. 
( Hi) for n large, 

e^ + rn{z) = {-IT {^-^)"''^\,n-,.{z)e^-' + ^ (i)) ' ^^'^^^ 
locally uniformly in C, where Cn is a constant that depends only on the interpolation 



points zj^"^ and such that 



c„ = 1 + C — , as n ^ oo. 



In particular, if pn = e{n)y/nwith e{n) which tends to as n tends to infinity, then {1.10) 
can be rewritten as 

e^ + rn{z) = i-ir [^y^^' W2n+i{z)e^'' (l + 0{e\n)) + 0(^)) . (1.11) 
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Remark 1.2 For the special case of bounded interpolation points, the error estimate ( |1.1 1[ ) 
agrees with the estimate (2.4) of [i2ril except for a minus sign that was incorrect there. 

Remark 1.3 For the theorem to be true, an assumption on the growth of p„ is mandatory. 
Indeed, if we allow p„ to grow linearly in the degree, that is p„ < cn with c > some constant, 
the theorem can be false. For instance, if c = 2tt, it suffices to consider the constant function 
rn{z) = 1 which interpolates the exponential at the points {±2i7rj, j = 0, . . . ,n} and does 
not converge to it as n tends to infinity. The particular case of shifted Pade approximants also 
shows that the theorem is false for linear growth, even for constant c smaller than 27r. Indeed, 
denote by Cq = 0.66274... the positive real root of the equation 

Vz^ + l + \og- ^== = 0. (1.12) 

1 + + 1 



Then, it follows from results in [1221 that shifted Pade approximants of degree n, interpolating 
at the point nc, where c is any real number with |c| > Cq, does not converge to e^. Still, we 
conjecture that Theorem 1 1 . 1 1 remains true if p„ < cn with c < cq. 

The next theorem describes the limit distributions, as tt, — )■ oo, of the zeros of the scaled 
polynomials P„ and Qn defined by 

Pn{,z) = pn{2nz), Qn{z) = qn{2nz). (1.13) 

For that, we need to introduce critical trajectories of the quadratic differential {z'^ + l)z~'^dz'^, 
defined by the condition 

Re r'^)^^.^0. ,1.,4) 



In (1.14) we assume that the square root has a branch cut along the path of integration and 
behaves like z at infinity. By (V + 1)+ we denote the + boundary value of the square root on 
that path of integration. An explicit integration of the differential form in the integral actually 



shows that condition (1.14) can be rewritten in the equivalent form Re {ri{z)) = 0, with ri{z) 



the expression in the right-hand side of ( 1 . 12 ). 



From the discussion in [[2T|. it follows that there are four critical trajectories, see Figure [T| 
We define 71 to be the critical trajectory connecting i with —i in the left half of the complex 
plane; the other critical trajectories are the mirror image 72 of 71 with respect to the imaginary 
axis, and the vertical half-lines (±i,±ioo). These curves determine three domains that we 
denote by Di^^o, Dq and -D2,oo as in Figure[Tj 

Next, we define two positive measures respectively supported on the curves 71 and 72, 
namely 

in s in s 

and, for a polynomial p of degree n, we denote by Up the normalized zero counting measure 

p = - 

n ^-^ 

p(^)=0 

Then, the following theorem holds true. 
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Figure 1: Critical trajectories satisfying (11. 141) and the domains Dq, Di^^o, -^2,00 



Theorem 1.4 As n tends to infinity, we have 

J^Pn^^P, J^Qn^f^Q, (1-16) 

where the convergence is in the sense of weak-* convergence of measures. 

The structure of the paper is as follows. In Section 2, we display a few basic properties of 
the rational interpolants and characterize them in terms of the solution of a specific matrix 
Riemann-Hilbert (RH) problem. In Section 3, we introduce different functions that are useful 
for the steepest descent analysis of the RH problem, which we perform in Section 4. From this 
analysis we derive in Section 5 our main results. Finally, in Section 6, we present numerical 
experiments in the case of interpolation points of order n. 



2 Rational interpolants and a Riemann-Hilbert problem 



As said before, polynomials p„ and g„ such that ( 1 .4 )-( 1 .5 1 hold true always exist. About 
uniqueness, we have the following simple proposition. 

Proposition 2.1 The irreducible form of the rational function Pn/^ln, associated to any pair of 



polynomials {pn, Qn) satisfying (1.4)-{1.5), is unique. 



Proof. Consider two pairs g„) and (p„, g„) satisfying ( 1 .4 )-( 1 .5 1. By taking away possible 



common factors, we get two pairs, each with coprime polynomials q'^) and (p^^, g^), 

max(degp^,degg^) = n - max(degpv„ deg g'^)) = n - da 
and sets Ji , J2 such that 

p'„iz) + q'^{z)e' =0(n,eJi('2-^r^))' ^1 C {0, . . . , 2n}, \J^\ > 2n + 1 - d,, 



(2n)- 



J2 C {0,...,2n}, IJ2I > 2n + l-rf2. 



Since (resp. g^) does not vanish at zf"^^, j e Ji (resp. Zj^"'\ j G J2), we may divide the 
above relations by q'^ and respectively, and deduce that 

\jeJlnJ2 / 

Since the degree of the polynomial in the left-hand side is at most 2n — rfi — ^2, we may conclude 
that Pn/Qn = P'n/On^ wMch implics uniqueness of the irreducible form of the rational function 
Pn/?n as asserted. □ 

Throughout, we will use the scaled interpolation points 

(2n) 

■ 2n , J u,...,/n. 



Note that, by assumption (1.8), we have 



|f"^l<c^, J=0,...,2n. 
Our main object of study will be the following Riemann-Hilbert problem. 

RH problem for Y 

Find a 2 X 2 matrix-valued function Y = F^^n) . £\^y^ ^ C^^^, with r„ a counterclockwise 
oriented closed curve surrounding the scaled interpolation points Zj^^\ j = 0, . . . , 2n, such that 

(a) Y is analytic in C \ r.„. 



(b) Y has continuous boundary values when approaching r„ from the inside (out- 

side) of r„, and they are related by the multiplicative jump condition 

Y4z) = Y4z) ' ' , (2.1) 



1 



with 



1 

V^{z) = 2z+-T^og{z-zf''^), (2.2) 

j=0 

(c) we have 

Y{z)z-''''' = I + 0{z-^), as 2 ^00, (2.3) 
where ^3 = f q '^i ) third Pauli matrix. 
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A motivation for the choice (|2.2[) of the potential Vn comes from the fact that 



-2nz 



z'PJz] 



where we set 

PJz) = 



-dz = 0, j = 0, . . . , n — 1, 



(2.4) 



VtJz) 



2n 

n 

3=0 



. - f ■') . 



These relations easily follow from ( |1.4 )-( [T3 ) and Cauchy formula. They can be interpreted as 
orthogonality relations for the polynomial Pn{z) on the contour r„ with respect to the varying 
weight e""^"^''). 

Next, we prove a proposition concerning existence and uniqueness of a solution to the RH 
problem, and relate this solution to the scaled polynomials Pn,Qn defined in (1.13), and the 
scaled remainder term 



E^{z) = Pn{z)e-^' + Qn{z)e^' = 0{n^{z)). 



Proposition 2.2 The following assertions hold true: 

( i) There is at most one solution to the RH problem for Y. For n large enough, a solution Y 
exists. 

(ii) Let n G N and Zj^"^ G Cfor j = 0, . . . , 2n. If the RH problem for Y has a solution, and if 
we write 



Pn{2nz) := Yii{z), 
qn{2nz) := a„(z)ri2(z) 



for z eC\ r„, 

for z outside ofTn, 



(2.5) 
(2.6) 



then PniQu are polynomials with degp„ = n, degg„ < n, and they satisfy the interpolation 
conditions ([73]). 

(Hi) Assume g„) is a pair satisfying (1.4)-{1.5) with 

degpn = n, degg„<n, 



and {pn-i, Qn+i) is a pair satisfying {1.1)-{1.2) with 

degpn-i <n-l, degg„+i = ra + l. 

Assume also that the normalizations ofpn and qn+i are chosen so that Pn{z) = pni^nz) and 
g„_i_i(2nz) are monic polynomials. Then, 



Pn{z) n-\z)Qn{z) 

yPn-i{2nz) Q:^^{z)q„,+i{2nz)^ 

Pn{z) n-\z)e-^'En{z) 

yPn-i{2nz) 1];^^ (2:)e~"^e„_i,„+i (2nz) 
solves the RH problem. 



Y(z) 



Y(z) 



z outside Tr 



z inside Tr 



(2.7) 



(2.8) 
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Proof. Uniqueness of a solution Y to the RH problem follows easily from Liouville's theorem, 
implying first that det Y{z) = 1 everywhere in the complex plane, and second, that the product 
YY~^(z), where Y is another solution, can only equal /, the identity matrix, since YY~^{z) 
has no jump and behaves like / + 0{l/z) at infinity. The fact that a solution exists for n large 
is a consequence of the steepest descent analysis to be done in Section |4j 

'1 



We now show assertion (ii). The first column of the jump matrix in (|2.1|) is ( ^ ) , so Yn 



and Y21 have no jump across r„, and they are entire functions. The asymptotic condition (c) of 
the RH problem tells us that Yn is a monic polynomial of degree n, which we denote by Pn, 
and that F21 is a polynomial of degree at most n — 1, which we denote by Pn-i- From the jump 



relation (2.1 1, it follows that 



27ri 



Pn(s)e 



-nV„{s) 



-ds. 



Calculating the integrals using residue arguments and the precise form ( |2.2[ ) of Vn, we find, for 
z outside r„, that Y12 is of the form 

Y^2iz) = n-\z)Q^iz) 



where Q„ is a polynomial which has to be of degree at most n because of ( |2.3| ). Similarly, we 
find that F22 is of the form 



Y22{z) = n-\z)q^^,{z), 



with g„_|_i a monic polynomial of degree n + 1. From the jump relation (2.1 ), it then follows 
that 



Y22(Z) 



-nVn{z) 



{pn-liz) + e^'^'qn+iiz)) , 



(2.9) 
(2.10) 



for z inside r„. Since Y12 is analytic inside r„, it follows from (|2.9[) and the definition (2.2 1 of K 



that the polynomials Pniz) = -Pn(^) and Qniz) = Qniw-) satisfy the interpolation conditions 



(1.5), as asserted. 



Assertion (iii) is easily checked. We leave the details to the reader. 



□ 



Corollary 2.3 For n large enough, polynomials pn and qn satisfying ( |i.4p -( [73| ) exist, are unique 
up to a normalization constant, and the interpolation problem is normal in the sense that the 
polynomials Pn and qn have full degrees. 



degpn 



n, 



degg„ 



n. 



Proof. For n large, we know from assertion (i) of Proposition 2.2 that the RH problem admits 



a solution Y which is unique. By assertion (ii) of the same proposition, the polynomials pn and 
g„ defined by (2.5) and (2.6) satisfy the interpolation conditions (1.4)-(1.5) with degp„ = n. 



Assume there exists another pair not a scalar multiple of the pair (p„, satisfying 

d, < d < n, 
rpn. 



( |1.4| )-( [T3] ). Then, by Proposition |2.1[ there exists a polynomial r, degr 
Pn = rpn, qn = ^qn- By Considering any other polynomial r of degree d, we get a pair 
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Qn = rqn satisfying (1.4)-(1.5), with degp„ = n, different from the pair {pn,qn)- Moreover, 
without loss of generality, we may assume that pn{2nz) is monic. Then, the matrix Y defined 
by 

^Pn{2nz) Q~^{z)qn{2nz 
Y(z) = \ 1 , z outside T 

Y2i{z) Y22{Z) 



ni 



IPn{2nz) . 
Y{z) = , z inside T 




is different from Y and, using assertion (iii) of Proposition 1X21 we see that it also solves the RH 



problem. This is a contradiction with the uniqueness of a solution Y . We may then conclude 
that, for n large, there exists, up to a multiplicative constant, a unique pair g„) satisfying 



(1.4)-(1.5), and also that degpn = n. The fact that, for n large, degg^ = ^ is a consequence 



of th e symmetry of our interpolation problem, namely that the pair {pn{z), qn{,z)) solves (1.4)- 
(1.5) w ith r e spec t to the interpolation points {-zp"^} if and only if the pair {qn{—z)^pn{—z)) 



solves (|l.4|)-(|l.5|) with respect to the interpolation points {—zf^^}. □ 



All the necessary ingredients to prove our convergence results, namely Theorem |1.1| and 



Theorem 1.4 are contained in the RH problem for Y . We will perform a rigorous asymptotic 
analysis of the RH problem for Y using the Deift/Zhou nonlinear steepest descent method [[H 

mum. 

Using this method, we will obtain existence and precise large n asymptotics for the matrix 
Y{z) everywhere in the complex plane. This allows us to obtain uniform asymptotics for the 
polynomials Pn and from which asymptotics for the original polynomials p„ and follow. 
Such asymptotics were obtained in [[2T]| for interpolation points satisfying ( 1.7[ ), and can be 
obtained similarly in our more general situation. 

3 Construction of the ^-function 

Assume that, for each n > 0, a suitable oriented curve 71 „, with endpoints a„ and 6„, is given. 
We want to construct a (7-function gn{z) satisfying the conditions 

(a) e^" : C \ 7i,„ C is analytic, 

(b) there exists G C such that 

^ 2n 

gnA^) + 9nA^) - 2;^ - -Y^^ogiz - zj'"^) + 2C = 0, for z e 7i,n, (3.1) 

where a branch of the logarithm is chosen which is analytic on 71 

(c) gn{z) = logz + (9(1) as 2; — 00. 
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This function will play an essential role in the asymptotic analysis of the RH problem for Y: it 
will enable us to transform the RH problem for F to a RH problem for T which is normalized 
at infinity (i.e. T{z) I as z ^ oo) and which has jump matrices which are convenient for 
asymptotic analysis. 

One can construct such a function gn with properties (a)-(c) for any smooth curve 71,^. 



Indeed, condition (3.1 ) is equivalent to 

2n 



= 2+^f] K^, (3.2) 



and the function 



with 

Rn{z) = {{z - an){z - bn)Y^\ Z G C \ 7i,n, (3.4) 

satisfies this condition together with the asymptotic condition g'^iz) ~ Mf the branch of _R„ is 
chosen which is analytic off 7i,„ and which behaves like z as z ^ oo. After a straightforward 
integral calculation, we get 

with 

1 ^ i?„(^'")) 

hn{z) = an + bn-2z + . (3.6) 

The (yf-function is then the multi-valued function 

gn{z) = / log(2; - s)i!n{s)ds, (3.7) 

where the branch cut of the logarithm follows 71 ,1 along (z, hn) and then it goes further to 
infinity. 

For the particular case of Pade approximants where all interpolation points z'f^\ j = 
0, . . . , 2n, are equal to 0, the above formulas become independent of n and reduce to 

1 z^-(^)z + R{0) 
g\z) = 1 + ^, (3.8) 

z zR[z) 



and 



= .1 , , (a + b-2z- , (3.9) 



27TiRj^{z) \ z / ' 
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with 



R{z) = {{z-a){z-b)y/\ 

and a and b two points independent of n. Now, recall the curve 71 that was defined before 
Theorem 1 .4 In the Pade case, we will choose that particular curve 71 in the definition of the 
(jf-function. Its endpoints i and —i satisfy ip(i) = ipi^—i) = 0, and we have 



RJz) 



TXIZ 



so that 



9{.z) = — [ log{z - s 
7r« J^, 



-ds 



(3.10) 



(3.11) 



becomes a complex logarithmic potential associated to a real measure, recall ( 1.14[ ). Actually, 
it is a positive measure by (1.15). Next, if we take a suitable curve 72 connecting +i with —i, 
lying to the right of 72, we have the important inequality 

Re(^+(z)+^_(2)-2(z + logz) + 2£) <0, for 72, (3.12) 

which is strict for z G 72 \ {±^}, see [i21i, Lemma 2.9]. We denote by V the closed contour 

r = 7iU72, (3.13) 

oriented counterclockwise. We return to the general case of complex interpolation points. From 
( |3.6| ), we obtain 

hja„) := 6n - On + - > , , =, (3.14) 



n 



and 



hn{bn) = an-bn + 



^ 2n 



n 



(3.15) 



In an ideal situation, we would choose a„ and 6„ in such a way that hn{an) = hn{bn) = 0, like 
in the Pade case. However, for a general set of interpolation points Zj^""^ bounded by ( |l.8 ), it 
is sufficient if hn{an) and hn{b„) are small enough. More precisely, we consider the following 
rescaling of the interpolation points. 



,{2n) 



2pn 



C(2n) 



t 



Pn ^ _^ 



n 



where p„ denotes the maximum of the moduli of the Zj^'^\ recall ( I.61. Hence the points Zj^"^ 
remain bounded with n, and we construct a„ and 6„ of the form 



(3.16) 
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in such a way that 



hn{an) = 0{t 



fc+l^ 



hn{hn) = 0{t 



A:+l^ 



as t — > 0, 



(3.17) 



with k a sufficiently large integer. We will see in Section 4.5 that ^ > 2^ is sufficient. Substi- 
tuting (3.16) in (3.14 )-( 3. 15 1, and expanding up to 0{t'^) gives us a system of 2k equations: the 
linear terms in t yield 



2ra 

^ ~(2n) 



. 2n 

I ~(2n) 



(3.18) 



and in general the (9(t'")-term in (3.14) (resp. (3.15)) gives an expression for am (resp. jSm) in 
terms of z[^"'^ for j = 1, . . . ,2n and in terms of aj, Pj for j = 1, . . . , m — 1. In particular, 



our system of 2k equations in the unknowns ai, . . . , and (3i, . . . , (3^ has always a unique 
solution. We can choose k in (3.17) as large as we want, but, as said above, k > ^ will be 
sufficient for us. 

Proposition 3.1 The points an and tend to i and —i respectively, as n tends to infinity. 



Proof. In view of (3.16), it is sufficient to prove that the coefficients ai and Pi, I = 1, . . . ,k, 
remain bounded as n tends to infinity. From (3.18), this is true when / = 1. 

For indices / > 2, we prove the assertion by induction. Let = —ian and bn = ibn, then 



hn{an) = -iF{an, bn), hn{bn) = -iF{-bn, -an), 



with 



^ 2n 

F{x,y) =x + y 

n ^-^ 



^/(a; + itzj){y — itzj) 
X + ifzj 



where we have dropped the superscript (2n) in the ^/s for simplicity. For / = 2, . . . , fc, it is 

easily checked that the coefficients of t' in the expansions of F(a„, bn) and F{—bn, —an) can 
be written as 



2n 2n 



2ai-2co- — Zj -Y] ij-, 



(3.19) 



2n 



n ^ 



, 2n 

n ^-^ ■' 



(3.20) 



where the Ci,. . . ,q and the c^i,. . . ,di are polynomial expressions in the ol\,. . . , . . ,Pi-\. 

Hence, the vanishing of (3.19)-(3.20) and the boundedness of the ij, j = 1, . . . , 2n, show 
inductively that the coefficients and / = 1, . . . , fc, of a„ and 6„ remain bounded as n tends 
to infinity. □ 
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Remark 3.2 The question arises whether the expansions (3.16) are convergent as — )■ cxd for 
small t. If this is true, we have 



hni lim a„) 

fc— ^-oo 



hni liin 



0. 



This would enable us to relax the condition ( 1.8 1 to 



Pn < CiTl, 

for a sufficiently small constant ci > 0. 



(3.21) 



(3.22) 



We still have to choose two curves 7i^„ and 72,„, each connecting a„ with 6„. These two curves 
will build up the closed contour r„ = 7i,nU72,n which appears in the Riemann-Hilbert problem. 
We take th e con tour r„ as a local deformation around the points i and —i of the contour F 
defined in (3.13). For that, we fix two sufficiently small disks U^"^^ surrounding ±i. From 
Proposition 



3.1 



for n sufficiently large, a„ and 6„ will lie in those disks. We let 71 „ and ^2,n 
respectively coincide with 71 and 72 for z outside U^'^\ and near ±i we can extend 7i^„ and 72,n 
arbitrarily, as long as they do not intersect and have no self -intersections. 



Finally, define 



^ 2n 

-2gn{z) + 2z+ -J^log 



for z outside of 71 „. For z G 71 it follows from (3.1 1 that 



(3.23) 



-9n,+ {^) + 9n-{z) 



and, together with ( |3.5[ ), this implies that 

-ds. 



hJs] 



Rnis) 



(3.24) 



wit h Rn and hn the functions respectively defined by (p4b and (3.6). The path of integration 



in (|3.24|) is in C \ (7i^„ U (U{^^"^})) and does not wind around any of the points 'zf'^' . The 



function (pniz) has logarithmic singularities at 2^^"^ j = 1, . . . ,n, and a branch cut starting 
at a„ which goes along 71 „ and then further to infinity. In the case of Fade interpolants, the 
function 0„ simplifies to 



ds. 



Proposition 3.3 Let U be a given neighborhood ofO. Then, we have 

(pn ^ (p, as n ^ 00, 

locally uniformly m C \ (71 U f/). Moreover, for any constant C < 0, there exists a neighbour- 
hood U ofO such that Re (0„) < C in U for n large enough. 

Proof. From the dominated convergence theorem, the facts that a„ tends to i and the points 
Sj^"^ tend to as n — )• 00 follows that 0„ tends to point- wise in C \ (71 VJU). By boundedness 
of the 0„ outside U, we derive that the convergence is locally uniform. The fact that Re (0„) 
has logarithmic singularities at the zj^"'' which tend to implies the second assertion. □ 
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Concerning the function 0, the following lemma about the sign of its real part will be useful in 
the sequel. 

Lemma 3.4 ( [|2T1 Lemma 2.9]) Let Dq, -Doo,i ^f^d -Dqo 2 be the open domains delimited by the 
four critical trajectories of the Pade case, as depicted in Figure [7] Then, the real part of cj) is 
negative in Di^^o U Dq, and it is positive in D2,oo- 



Note that, by ( 1.14l, the real part of vanishes on 71, 72 and the vertical half-lines (±2, ±ioo). 



4 Steepest descent analysis of the RH problem 

As usual, the steepest descent analysis consists of a number of transformations. 



4.1 First transformation Y ^ T 



Define 



(4.1) 



Note that r„, which is fixed outside U^^\ will surround the scaled interpolation points " for 
n large by ( |1.8[ ). 

RH problem for T 

(a) T : C \ r„ ^ is analytic, with r„ = 7i,„ U 72,n, 

(b) T has continuous boundary values when approaching r„, and they are related by 

T+{z)=T4z)Jt{z), 

with 



Mz) 



~n{gn,+ {z)-gn,-{z)) ra(g„,+ (z)- Vn (z)+2£„) 







Dn{gn,+ (z)-g„-{z)) 



(c) T{z) =1 + 0{z-^) as ^ ^ 00. 



Making use of the function 0„ defined in (3.23), we obtain 



M^) = { 



gn</>„,+ (z) giy„{2) 

g gn</>„,_(z) 
'1 QWn{z)^-n4'n(z)' 

1 



with 



for z G 7i,n, 

for Z G 72,n, 



WJz) = -loefz - z 



(4.2) 



14 



4.2 Opening of the lens S 

On 7i „, the jump matrix for T can be factorized: 



e^"(^) 

gn0„,_(2;)g-W„(2) I I g-VK„(z) g I 1 gn0„,+(2;)g-W„(2) 



(4.3) 



The jump matrix Ji can be continued analytically to a region to the left of 7i,„, and J3 to a 
region to the right of 71,^ (simply by replacing 0„_+ and 0„ _ by This enables us to split 
the jump contour 71 „ into three distinct curves 7J „, 71 „, 7^'^ (from left to right) connecting +i 
with —i: we call this the opening of the lens. Define 

{T{z), for z outside the lens-shaped region, 

T{z)Ji{z), for z in the left part of the lens, (4.4) 
T{z)J3{z)^^ , for z in the right part of the lens. 

RH problem for S 

(a) 5 : C \ r„ ^ is analytic, with r„ = 7i,„ U 72,„ U -f^^ U 7;;„, 

(b) 5* has continuous boundary values when approaching r„, and they are related by 

S+{z) = S-{z)JTiz), for2;G72,n, 

S+iz) = S^iz)Ji{z), forzey^^^, 

S+{z) = S-{z)J2iz), for2;G7i,„, 

S^z) = S.iz)Jsiz), iorze^l^, 

(c) S{z) = 1 + 0{z-^) as 2 ^ cx). 

We choose 7J „, 7"^ to lie sufficiently close to 71 „, and in any case in the region where 0„ is 

analytic, away from the scaled interpolation points Zj^"'\ The jump matrices Ji, J3, and Jt 
decay on the contours 7^ „, 7"^, and 72,n respectively as n — 00, except in two small but fixed 



neighborhoods U^"^^ of ±i. This follows from Proposition 3.3 and Lemma 3.4 together with 
the fact that Wn{z) is uniformly bounded on the jump contour. The only jumps that survive the 
large n limit, are the jump on 71 „ and the jumps in the vicinity of ±i. In this perspective one 
is tempted to believe that the leading order asymptotic behavior of S, away from ±i, will be 
determined by the solution to a RH problem with jump J2 on the curve 71 „. Nevertheless, some 
substantial analysis remains to be done to turn this into a rigorous argument. 

4.3 Outside parametrix 

Consider the following RH problem, on the contour 71 „ connecting a„ with 6„. 
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RH problem for P(°°): 

(a) P(°°) : C \ 7i,„ ^ C^^2 is analytic, 

(b) Pt\z) = p[-\z) ' Q J , for ^ e 7i,n, 

(c) P(°^)(2) = J + C(z-i), as 2^00. 
A solution to this problem is given by 

p(°^\z) = ( ^^-^ ] NDn^zy^ (4.5) 



for z G C \ 7i^„, with 



z — ar. 



1 , 



iV = 4= ( \ ! ) e-^™^ (4.6) 



and Dn{z) is the Szego function which is analytic and non-zero in C \ 7i,„, and satisfies 

An explicit expression for Dn{z) is given by 

D„W=exp|'^/ „ 7"';> (4.7) 

27rz J^^^^ Pn,+ (s)(s - 2;) 



with i?„ defined by (3.4). This outside parametrix p(°^) = P^°°\z;n) will determine the 
leading order asymptotics of S{z) for 2; in C \ f/^^\ as n — t- 00. In U^'^\ we need to construct 
local parametrices P*-^^ that determine the leading order asymptotics of S in those disks. This 
is the goal of the next section. 



4.4 Local Airy parametrices 

We will construct local parametrices in the regions U^"^^ surrounding ±i. These parametrices 
P = P^^^ should be analytic in U^^^ \ (7^ „ U 7i^„ U 7"^ U 'y2,n), see Figure |2| and they should 
have exactly the same jumps as S has on fl (7^ „ U 7i^„ U 7"^ U 72,ri)- In addition, we 
aim to construct the parametrices in such a way that P^'^^z) P^°°\zy^ is as close as possible 
to the identity matrix on dU^'^\ As is common for the construction of local parametrices near 
the points where the lens closes, we will build P using the Airy function. If our function 
(pniz) would behave like c{z ± a^)'^/^ as 2; ^ a„ and as c(z ± &„)^/^ as 2; — )• bn, this would be a 
standard construction as in [l5l|6l|71. Unfortunately, this is only the case if /i„(a„) = hn{bn) = 0. 
Therefore we will need some technical modifications to have suitable parametrices. The reason 
why we are still able to construct Airy parametrices, is that /i„(a„) = /in(^n) = 0{t'''^^) as 
n — > 00. 
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Figure 2: Contours in the neighborhood C/*^ ^ around the point —i 

4.4.1 Airy model RH problem 

Define 

% = % (C)=cu^Ai(a;^C), J = 0,1, 2, 
with u! — and Ai the Airy function. Let 



where yj denotes the derivative of yj with respect to (. Since the Airy function is entire, each 
Aj is an entire matrix function. Furthermore, using the identity yo + yi + y2 — 0, it follows that 



A(C)=^4(C) 



1 1 

1 



A(C) = ^(C)(J J) 



^2(0=^3(0 
^3(0 = ^4(0 



1 
-1 

1 

1 1 



From the asymptotics of the Airy function and its derivative, 

Ai(c) = ^r"vi<-(i+c(J^)), 



(4.8) 

(4.9) 

(4.10) 
(4.11) 

(4.12) 



Ai'(C) 



-1 
2^ 



l + O 



^3/2 



(4.13) 
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Figure 3: The four sectors (rj)j=i_,..^4, and the contours and jumps for the matrix A{Q. 



as C — oo with | arg^| < tt, follows that 

as C — > oo in the sector Tj defined by 

Ti = {CgC: <argC<7r}, 

f2 = {CeC:^<argC<y}, 

' — ' Stt tt 

r3 = {CGC:-y <argC<--}, 

r4 = {CeC:-7r<argC<^}. 



1,...,4, 



(4.14) 

(4.15) 
(4.16) 
(4.17) 
(4.18) 



Let Fj, j = 1, . . . , 4, be the sectors delimited by the four rays of argument 0, tc, as 
shown in Figure [3] Then, it follows from what precedes that the matrix A such that 

A{C) = MC), Cer„ j = i,...,4, 

admits the jumps shown in Figure |3] and has the asymptotic behavior 

A{C) = C^^N [I + O (C~'/')] e-^f'^'^^a, as C ^ oo, (4.19) 



with the constant matrix defined by (4.6 1 



4.5 Construction of the parametrix in U''' 



We search for functions /„, Sn in ^ such that the function (pn, defined by (3.23), can be 
expressed as 



Mz) = -Uz)'/^ + 2s„(z)/„(z)i/2 mod 2m, 



for z G U^'K 



In view of the integral expression (3.24) of (j)n{z), we therefore define /„ by 



hn{s) - K{K 



Rn{s) 



-ds. 



(4.20) 



(4.21) 
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and Sn{z) by 



hn{K 



-ds. 



Then /„ and s„ are both analytic functions in U^^\ and we have 



(4.22) 



fniK 



0, 



3/2 



\KXbn)\ 



A,•+l^ 



2a/|o„ - a„| 
as n — 7- oo, uniformly for z E U^~\ 



(4.23) 
(4.24) 



For n sufficiently large, fn + s„ is a conformal map from ^ onto a neighborhood of 0. In the 
Fade case, we have 



ds. 



and it suffices to consider the function / defined by 



4 



2/3 



which is analytic in a neighborhood of —i, and where the 2/3rd power is taken so that f(z) is 
real negative for 2 G 71. Then, the three contours 72, 7^ and 7" are chosen so that they are 
respectively mapped by / on the real positive semi-axis, and on the rays of argument — 27r/3 
and 27r/3. In our situation, and for n large, the map /„ + Sn does not send 71 „ exactly on the 
negative real line, see Figure|4j Choosing 71 „ as an arc which tends to 71 as n tends to infinity, 
and since /„ + s„ converges uniformly to / in f/^") as n tends to infinity, we get that 71 „ is 
mapped to an arc Ai „ which tends to the negative real line. Then, the contours 7^ 7"^ and ^2,n 

can also be chosen so that they are mapped by /„ + s„ to contours A'^ „, A" „ and A2,ri tending to 
the three other rays of the usual Airy parametrix. 

For E^^^ some matrix-valued analytic function in U^^^ that we will define in the sequel, let 



P^-\z) := E^-\z)A„{n'/'iUz) + s^iz)))e'^^-^^^^h 



(4.25) 



with 



1 4, 



and Sj^n is the image by the conformal map + s„ of the region i?j „ inU'^ \ as indicated in 
Figure|4j As a consequence of (4.8 )-(4.1 1 1 and the right multiplication with the two exponential 



factors in ( |4.25[ ), it follows easily that 

Pi-\z) = P[-\z).Mz), 
Pi^\z) = P[-\z)Mz), 



Pi~\z) 



P^-\Z)J2{Z), 



for Z e 72,n, 

for z G 7l_„, 
for z G 7i,n, 
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fn{z) + Sniz) 



Figure 4: Images of the domain ^ and the contours by the conformal map /„ + s„. Inside 
the image domain, contours of the usual Airy parametrix are drawn in dotted lines. 



Pi-\z) = Pr\z).Mz), for zeyi^. 

Now our main concern is the matching of P^~^ with p(°°) at dU^~\ Suppose that z is in 
region Rjn and on the boundary of U^~\ then one can verify that n^/^{fn{z) + Sn{z)) lies in 
region Vj defined in (4.15 )-(4.18 1, if n is sufficiently large. Consequently we can use (4.14) for 
z E dU^^\ and we obtain 

P^-\z) = E^-\z){n"\Uz) + .„(^)))-^iV 

X [/ + 0(n-i)]e-i'^(^"(^)+^"(^))''''^^et^"(^)'^^e-'^'^^ (4.26) 



Substituting (4.24), we find by (4.20) and t < cn'°' that (4.26) simplifies to 



0{n 



l-2(fc+l)a 



0-3 



e 2 

(4.27) 



For k > have 

P^-\z) = E'^-\z)n-'^ fn{z)-'^ N [I + C(n-2S)J e-^"^ (4.28) 
with a = min{a, 1/2}. Since we want p(~^ to match with the outside parametrix, we define 

E(-)(z) = p(°°)(2)e^"^Ar"V„(2)^n?, (4.29) 

and we need to check that E^^'^ is analytic in \ since, i f not, the jump relations for P^^^ 
would be violated. By (4.5 ), it is easily checked that, indeed, (4.29 ) defines E^^^ analytically in 



U 



(-) 



Since the function Wn defined by (4.2 ) is bounded on dU^'^\ we have 



P^-\z)P^ 



oo). 



I + 0{n 



-2a\ 



for z G dU'^ \n ^ oo. 



(4.30) 
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Figure 5: Contour S/?. 



A similar construction works for the local parametrix P*^+) near +i if we let P(+) be of the 
form 

PW(^) := EW(z)a3A„(-n2/3(/„(^)+s„(z)))a3e"^'^^e-'^'^^ for z in region P;.„, 

(4.31) 

with region R[ „ being the one outside the lens and to the left of it, and the regions R'2 „, -R3 „, i?4 , 
occur in order when turning around i in counterclockwise direction. One can mimic the con- 
struction of P*^^ with minor modifications such that 

p(+)(z)P^'^\z)-'^ = I + 0{n~'^^), for z e dU^+\n^ 00. (4.32) 



4.6 Final transformation 

Now we define 



R{z) 



S{z)P^°^\z)-\ for^GC\f/W, 



S{z)P^ 



(4.33) 



Then R is analytic in C \ with as shown in Figure [5} and it tends to / as 2; — 00. On 
S/j, we have 



R+{z) = R.{z) (/ + 0(n-2")) 



(4.34) 



unifor mly as n — )■ 00. The 0{n ^'^) error term is present for z E dU^^^ because of (4.30) and 
(4.32). OnSij\9f/(^) the ^umps are even exponentially small as n — )• 00, see the discussion 
about the RH problem for S at the end of Section |4.2[ Standard estimates in RH theory show 

as n — )■ 00, (4.35) 



that (|4.34[) implies existence of the RH solution R for large n, and the asymptotics 



R{z) =I + 0{n 



-2a\ 
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uniformly for 2; G C \ Sr. Reversing the explicit transformations (4.33 1, (4.4), and (4.1 ), we 
have existence of Y for n sufficiently large, and we can also find asymptotics for F as n — > 00. 
For z outside the lens-shaped region and outside U^'^\ we obtain 

= ((1 + 0{n-^^))P^^\z) + 0{n-^")PiT\z)) e"'"^'\ (4-36) 

and 

= ((1 + 0{n-^^))Pif{z) + 0{n-^^)Pi^\z)'^ ^-n9„(z) ^-2ne„ ^ (437) 



5 Proof of the main results 

Let 



2 \ \ z — i J \z — i 



1 ( f z — b„\^^^ f z — br,^ 



" 2 \\z — an J \z — a 

where the 1/4-roots are defined outside of 71 and 71 „ respectively, and tend to 1 as 2; tends to 
infinity. Note that, since a„ = i(l + 0(n~°')) and = — «(1 + 0(?t.~")), as tt. — )• 00, we have 

A^\z) = r^^\z){l + 0{n-'')), n^oo, (5.1) 

locally uniformly in C \ 71. Also, since 2^^"^ = 0{n^°'), we have that 

= D(^)(l + 0(n-")), n^oo, (5.2) 

locally uniformly in C \ 71, where 



2m J^^ R+{s){s - z) 

We first prove the following proposition which gives asymptotic estimates for the polyno- 
mials Pniz), Quiz) and the error function En in the complex plane, respectively outside of the 
curves 71, 72 and {±i, ±00). 

Proposition 5.1 As n —i- 00, we have 

p^(z) = rW(2)D-i(z)e"^"(^)(l + 0(n-°)), (5.3) 

uniformly for z in compact subsets ofC\ 71, 

' -r^+\z)D-\z)e''^3n{z)-2z)(^l ^ 0(n-°)), 2 G L'o U 71 \ {i, -i}, 

Quiz) = { _ (5.4) 

-ir]„(z)r(-)(2)L)(z)e-"(S"(^)+2'")(l + 0(r2-°)), zeC\Do, 
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uniformly for z in compact subsets ofC\ 72. Furthermore, we have 

'r(+)(2)D-i(2)e"(f"(^)-^)(l + 0(n-")), z E ^1,00 U 71 \ {i, -i}, 



En{z) 



(5.5) 



uniformly for z in compact subsets ofC \ {±i, ±ioo). 



Proof. The outside parametrix P^°°\z) dep ends on the endpoints a„ and 6„. We know that 
Qn —^i,bn—>—iasn—^ 00. In view of (4.5 1, we can then conclude that the entries of p(°°) (z) 



have modulus uniformly bounded below and above in compact subsets of C \ 71, as n — > 00. 



Since Pn{z) = Yii{z), we thus get from (4.36) that, locally uniformly, 



Moreover, it follows from (4.5 1 that 



PiT\^) = rl^\^)D-\z), 



which implies (5.3) because of (5.1 ) and (5.2). 

For Qn, three different cases need to be considered. First, we assume z E C \ Dq. Then we 



can assume that the curve 72 is such that z lies outside the contour T. From (2.7 ) and (4.37 ), we 
get 

Qn{z) = Yu{z)Qn{z) = P^f{z)Q4z)e-''^-^^^e-''''-{l + (9(^-2-)). 



From the fact that 



Pifiz) 



-iri7Hz)Dr, 



and using ( |5.1[ )-( |5.2[ ), we obtain the second estimate in ( |5.4| ). For z E Dq, we use the fact that 



Qn{z) 



-""'EJz) 



-2nz 



PJz) 



(5.6) 



and we need to find out the dominant term as n gets large. From ( |5.3[ ), we have that, as n tends 
to infinity. 



1 



n 



log |e' 



-2nz 



Pniz) 



Re{gn{z)-2z) + 0{n 



-1^ 



Similarly, from the second formula in ( |5.5| ), which we will prove next and independently, we 
obtain 

- log |e-"^E„(z)| = Re {~gn{z) - 2/„ + - logfi„(z)) + 0{n-'). 
n n 

The difference between the right-hand sides of the two previous estimates equals —Re (0„(z)) + 



@{n ^) which, in view of Proposition 3.3 and Lemma 3.4 is positive, locally uniformly in 



Dq, for n large. This implies that the dominant contribution in (5.6) comes from the term 
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— e 



"^"^P„(2). Hence, the first estimate in (5.4) for z E Dq follows from (5.3). It remains to 
check that the previous estimate still holds true when 2; G 71 \ We do not give the 

details here and we simply refer to [|271 Theorem 2.10] where a proof of a similar assertion is 
given. 

For the error function En, and for 2; G -Do. we have 



En{z) 



12U 



P'lf{z)nn[z]e 



n{z-gn(z)) -2ntn 



which leads to the second estimate in ( |5.5| ) when z G -Do- When z is in Di^^o or -D2,oo. we use 
that 



EJz] 



Pn{z)e 



+ Qn{z)e' 



and find out which term is dominant in the sum. Using the previous estimates (5.3) and (5.4), 
along with Proposition |3.3| and Lemma 3.4 it turns out that P„(2;)e^"^ dominates in -Di 00 while 
Qniz)e"'^ dominates in -D2,oo- Then, (5.3) and (5.4) are used to derive (5.5) in -Di oo and -D2,oo- 
Finally, one checks that these estimates are also valid on the curves 71 and 72. □ 



5.1 Proof of Theorem 1.4 



Proof. Applying Rouche theorem on a circle sufficiently large to contain the curve 71, we 
obtain, in view of the asymptotic estimate ( 5.3[ ), that for n sufficiently large, the difference 
between the numbers of poles and zeros of P„ outside of the circle equals the corresponding 



difference for the product of functions in the right-hand side of (5.3 ). This product has no zero 



(note that, in view of ( 3.7 ), Re (gn) (z) is lower bounded if z stays at some distance from 71) but 
e"'^"-^^^ has a pole of order n at infinity, since gniz) = log z + 0(1) as 2; — > 00. Since P„ has a 
pole of multiplicity n at infinity, we may thus conclude that P„ has no zero outside of the circle 
for n large enough. A similar argument using Rouche's theorem shows that, for any compact in 
C \ 7i, there exist tiq such that P„ has no zeros in that compact for n > uq. 

Now, let us consider the sequence of counting measures z/p„, n > 0. We already know that, 
for n large enough, all these measures are supported inside a fixed compact set. From Helly's 
selection theorem, we may thus select a subsequence of up^^ converging in weak-* sense to a 



measure u. Besides, from ( |5.3[ ), it follows that, as n — )• 00, 

1 



n 



\og\Pniz)\ = Reigniz)) + 0in- 



zeC\^u 



and from the dominated convergence theorem, we see that the sequence of functions gn{z) tends 
to g{z), defined in (3.1 1 ), point-wise in C \ 71. Hence, as n — )• 00, 

1 



n 



log \Pn{z)\^ Re {g{z)), zeC\^i, 



or equivalently. 



log \z - s\dup^{s) Re {g{z)), z G C \ 71. 

Since up^ tends to u and P„ has no zero inside any compact set of C \ 71 for n large enough, 
the above integral also tends to J log \ z — s\diy{s), point- wise in C \ 71, as n — )■ 00, so that 



log \z — s\di'{s) 



log \z - s\dfip{s), zeC\-fi, 
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where we recall that the function g is the complex logarithmic potential associated to the mea- 
sure djip. Since 71 has two-dimensional Lebesgue measure 0, the unicity theorem [13, Theorem 
11.2.1] applies, showing that v and fxp are equal. Since yUp is the only possible limit of a weakly 
convergent subsequence, the full sequence up^ converges weakly to iip. 

The similar result for the sequence of counting measures z/q„, n > 0, follows from the 
symmetry of our interpolation problem, already mentioned at the end of the proof of Corollary 
[Q □ 



5.2 Proof of Theorem [LT 

Before we start the proof, we need a few preliminary results. 
Lemma 5.2 Let -u G C \ 7„,i. Then, we have 
Rjz) f log{u - s 
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7n,l 



Rn,+ {S){S - Z 



1 ( {Wi{u) - W2{z)){Wi{z) - W2{u)) 

-as = - log 



(5.7) 



1 
27ri 



7n.,l 



{an + hn- 2s) \og{u - s) 

Rn,+ {s) 



ds 



(an - hn) 



2 2w2(z) + an + hn 



where 



Wi{s) = -S + Rn{s), W2{s) = -S - Rn{s), S G C \ 7i,„. 



(5.8) 



(5.9) 



Proof. We will not give the details for the proof of these formulas. We just mention that the 
computations can be performed by using the change of variables s ^ w, where s and w satisfy 
the algebraic equation 

w"^ + 2sw + {an + bn)s - a„6„ = {w - wi{s)){w - W2{s)) = 0, 



and then applying the Cauchy formula in the w-plane. 



□ 



Proposition 5.3 The functions Qn, and the constant 2£„ admit the following explicit expres- 
sions, 



1 



(a„ — br. 



9n{Z) 



Di{z) 



2 2w2{z) + a„ + hr. 



2n 



1 

2n 



{W2{'zf"^) - Wi{z)){wi{zf"^) - W2{z)) 



2w2{'zf''^) + an + b„ 



-2w2{z) -an-K 



2L 



{w2{z^q"'^) - Wi{z)){wi{z'q''^) - W2{z)y 

^ I f 2w,{f''^) + an + bn 

j=l \^W2{Z) ')+an + bn 



(5.10) 
(5.11) 

(5.12) 



Proof. The proofs of (5.10) and ( |5.1 1[ ) simply follow from the previous lemma together with 
the expression (3.5)-(3.7) of gn and the expression (4.7) of Dn- The proof of (5.12) follows by 
plugging (|5.10|) into (3.1) and performing a few calculations that we do not detail. □ 
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Proof of Theorem Assertion (i) is a consequence of the strong asymptotics obtained in 
Proposition |5.1 1 In order to prove (1.9) we need to evaluate the asymptotic estimates (5.3 1 and 
( [54 ) with z replaced by 2;/2n where z is a fixed complex number. We get 



Qn 



Z 

2n 

z 

2n 



rW(0)D-^(0)e"«"(^)(l + 0(n-")), 
-r(+)(0)D-^(0)e"^"(^)e-"(l + 0(n-")). 



(5.13) 



(5.14) 



Since the normalization chosen in Theorem |1.1| is such that g„(0) = 1, we can deduce from 
(|5T3l)-(|5T4l) that 



p^(z) = -e"(5"(5^)-5"(°))(l + 0(n-")), 
As gn tends to g uniformly in a neighbourhood of 0, we have that 



ngn 



z 

2n 



ngn{0)+g'M^ + & 



(5.15) 
(5.16) 

(5.17) 



where the 0(l/r2) term is uniform in z. Finally, g'^iO) tends to 
9(0) - 



1 f iV^^] 



ITT 



-ds. 



71 



which, by Cauchy theorem, is easily computed to be 1. Hence, by plugging (5.17 1 into (5.15 )- 
(5.16), the limits ( 1 .9 ) follows, and the convergence is locally uniform in C. 
For the third assertion, note that 



P 

e^ + r„(z) = e"+ "^^n. 



2n' 



Qn( 



Qni 



2n' 



Substituting ( |54| ) and ( [53] ), we obtain by ( [5lJ ) and 

r(-\0) W2n+l{z) 



e + r„{z) = I 



r(+)(0) {2n] 



2n+l 



D(0)2e2"e-2"^"(^)e-2"^" (l + 0(n-°)) . (5.18) 



Evaluating r'^^\Qi) and -D(O) using Proposition 5.3, and expanding ngn{^), we find 



e + r. 



(5.19) 



where we have also used that g'{0) = 1. From (5.10) and (5.12) we can derive an explicit 
expression for 2gn{0) + 2in, namely. 



2n 



2^„(0)+2C = -2^/^+— J2 log 
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(5.20) 



Moreover, from ( |3.16[ ), ( |3.18[ ) and the fact that the coefficients aj and /3j are bounded with 
respect to n follows that 



2n 



k=l 



Also, recalling the definitions (5.9 1 of Wi and W2, we deduce that 

2n 



A:=l 
2n 



k=l 



Plugging these estimates into ( |5.20p we get 

2(?„(0) + 24 = -2 + log(-4) + Oit") 



which together with (5.19) shows (1.10) and the assertion about the constant c^. 



□ 



6 Numerical experiments with interpolation points of modu- 
lus as large as the degree of the interpolant 

In this section, we are interested in the location of zeros and poles of rational interpolants 
associated to interpolation points whose moduli are comparable to the degree of the interpolant. 
For the particular case of shifted Pade approximants pn"^ /qn "'^ interpolating the exponential 
function at we have the simple relation 



pt\z) 
qt\z) 



qn\z - in 



where pn"* {z) / qli'' (z) denotes the usual Pade approximant interpolating at 0. Hence, we know 
at once that the distributions of poles and zeros follow exactly the shift of the interpolation point 
in ■ In particular, the limit distributions of the scaled (by 2n) zeros and poles are modified if 
and only if in/n does not tend to 0, as n tends to infinity. 

In Figure [6} we consider the case of 2-point Pade approximants with two real symmetric 
interpolation points of equal multiplicities. As the points approach the zeros and poles of the 
usual Pade approximant (or equivalently, in the scaled situation, the critical curves 71 and 72), 
we can see how the distributions of zeros and poles are modified. Clearly, a repulsion takes place 
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-100 -50 50 100 -100 -50 50 100 



Figure 6: Zeros and poles (circles in the left and right-half planes respectively) of 2-point Pade 
approximants of degree (51,50) with two real interpolation points (diamonds) of multiplicity 
51 located at {—50, 50}, {—65, 65}, {—85, 85}, { — 100, 100}. For comparison, the zeros and 
poles of the Pade approximant of degree 50 are shown with dots. 



between the interpolation points and the zeros and poles of the approximants. In Figure |7| we 
consider the case of interpolation points regularly distributed on a real segment. As the segment 
approaches the zeros and poles of the usual Pade approximant, we can again observe how the 
distributions of zeros and poles are modified. Finally, in Figure [8} the case of interpolation 
points regularly distributed on a circle is depicted. The zeros and poles of the interpolants seem 
to be pushed away as the circle intersects the limit distributions corresponding to the usual Pade 
approximant. 

From a theoretical point of view, polynomials whose roots are the above zeros and poles 
still satisfy the orthogonality relations ( |2.4 ). We note that, in the corresponding potential (2.2), 
the sum of log terms becomes preponderant as the interpolation points grow faster with n. 
This should account for the modification in the zeros and poles distributions of the rational 
interpolants, that we observe in our experiments. In any case, it would be interesting to study 
in more detail the interaction between the interpolation points and the zeros and poles of the 
approximants. 



28 



100- 




-100- 



100 



0- 



-50 



-100 



-100 


-50 





50 


100 




■o 

/ 










\^ 








-100 


-50 





50 


100 



Figure 7: Zeros and poles (circles in the left 
interpolants of degree 50 corresponding to 101 
real segments 60/, 72.5/, 87.5/, 110/, / = [-1 
Pade approximant of degree 50 are shown with 
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and right-half planes respectively) of rational 
interpolation points regularly distributed on the 
, 1] . For comparison, the zeros and poles of the 
dots. 
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